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The static Coulomb potential of Quantum Electrodynamics (QED) is calculated in the 
presence of a strong magnetic field in the lowest Landau level (LLL) approximation us- 
\ ing two different methods. First, the vacuum expectation value of the corresponding Wilson 

loop is calculated perturbatively in two different regimes of dynamical mass mdyn., 
|q^| <C TT^dyn ^ l^^l ^^'^ ^dyn Iq^I ^ where q|| is the longitudinal components of 
the momentum relative to the external magnetic field B. The result is then compared with 
^ , the static potential arising from Born approximation. Both results coincide. Although the 

QJ , arising potentials show different behavior in the aforementioned regimes, a novel dependence 

' on the angle 6 between the particle-antiparticle's axis and the direction of the magnetic field 



(N 



is observed. In the regime |qjj| ^ ™dyri ^ k^l, for strong enough magnetic field and de- 
^ i pending on the angle 9, a qualitative change occurs in the Coulomb-like potential; Whereas 

\l , for = 0, TT the potential is repulsive, it exhibits a minimum for angles 6 €]0, n[. 
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INTRODUCTION 



^ ' Chiral symmetry plays an important role in elementary particle physics. It has been known ilj for some 
H ' 

' time that QED, in addition to the familiar weak-coupling phase, may have a nonperturbative strong 



coupling phase, characterized by spontaneous chiral symmetry breaking [4]. The existence of this new 
phase was exploited in a novel interpretation of the multiple correlated and narrow peak structures 
in electron and positron spectra observed at GSI several years ago [3]. According to this scenario the 
electron-positron peaks are due to the decay of a bound electron-positron system formed in the new QED 
phase induced by a strong and rapidly varying electromagnetic field which is present in the neighborhood 
of a colliding heavy ions. It is therefore of great interest to investigate whether background fields, such 
as constant magnetic fields, can potentially induce chiral symmetry breaking in gauge theories and lead 
eventually to the formation of chiral symmetry breaking fermion condensate (i/jip) and a dynamically 
generated fermion mass. 

Indeed, the magnetic catalysis of dynamical chiral symmetry breaking has been established as a 
universal phenomenon in 2-|-l and 3-|-l dimensions 0, Q]- According to these results, even at the weakest 
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attractive interaction between fermions, a constant magnetic field leads to the generation of a fermion 
dynamical mass. The essence of this effect is the dimensional reduction D ^ D — 2 in the dynamics 
of fermion pairing in a magnetic field, which arises from the fact that the motion of charged particles 
is restricted in directions perpendicular to the magnetic field. It is believed that at weak coupling this 
dynamics is dominated by the lowest Landau level (LLL). The magnetic catalysis is not only interesting 
from purely fundamental point of view, but it has potential application in condensed matter physics [5| 
and cosmology 0]. 

In this paper we are interested on the static potential between the particle and antiparticle in the 
presence of a strong constant magnetic field. The potential produced by a point electric charge placed 
into a constant magnetic field is recently calculated in [7]. It is shown that the standard Coulomb law 
is modified by the vacuum polarization arising in the external magnetic field. Here, since the vacuum 
polarization component, taken in one-loop approximation, grows linearly with the magnetic field, a scaling 
regime occurs in the limit of infinite magnetic field. The scaling regime implies a short range character 
of interaction, expressed as Yukawa law 



, , ae"*^^^ l2a\eB\Nf , , 

y(x) = with \ ^ ^^-^^ 



where R = |x| and My is the photon mass. 

In the present work, we will determine anew the potential V{x) between the particle-antiparticle pair 
in the presence of a constant magnetic field in the LLL approximation. Although our method differs 
essentially from the analytic methods used in 0] leading to (jl.ip , our result is indeed consistent with this 
potential. 

To determine the fermion-antifermion potential, we will first compute perturbatively the vacuum 
expectation value (vev) of a Wilson loop of a static fermion-antifermion pair for large Euclidean time.^ 
In 0], the same perturbative method is used to determine the quark-antiquark potential of ordinary 
QCD. At one-loop level it is shown to be 

/ lle^ R\ 

where C = |tr (A^A'^) is the trace over the product of Gell-Mann matrices A"^, and is the UV cutoff 
parameter. Using the same idea, we will determine the potential of a point charged particle in the external 
magnetic field. Eventually, we will compare our result with the modified Coulomb potential from a semi- 
classical Born approximation. In the regime of LLL dominance, we will consider two different regions of 
dynamical mass, |qy| <C w-^y^ *C \e-B\ and m^^^ ^ |qy| <C |e-B|, separately. Here, qy is the longitudinal 
component of the momentum with respect to the direction of the external magnetic field. We will show 
that the potential in the regime |q^| <C ?ti-^j,„ ^ |e-B| has the general form of a modified Coulomb 



potential 



v,{R,e) = --lMa,e)-^^+ h -th 7 = ^^^, (1.3) 



dyn. 



^ Here, to determine thepotential in the LLL approximation, we will use the full photon propagator in the presence of a 
strong magnetic field [l^ in the LLL approximation. In [l^, it is shown that the full photon propagator depends on the 
dynamical mass mdyn., which is calculated nonperturbatively in [3]. Thus although our method is a perturbative one, the 
result depends automatically on a parameter which is determined nonperturbatively. 
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and in the regime 'rn'^y^ <^ |q|| <^ \eB\ has the form of a Yukawa- hke potential 

V2{R,e) = - ^^_^^^^^^^ , with M,exe)^g{e)^^-^, (i.4) 

where Ai,i = 1,2,3 in ()I.3P and g{9) in ()I.4p will be calculated exactly in Sec. III. In (jOj) as well as 
p.4p . ^ is the angle between the particle-antiparticle axis and the direction of the magnetic field. Up to 
this explicit novel dependence on the angle 9, the potential V2{R,9) from (|I.4p is comparable with the 
potential p.ip from [7]. As a consequence of this 6'-dependence, the effective photon mass M^q (9) in 
p.4p is, in contrast to the photon mass in (jl.ip . a function of 9. 

In the regime |q|| <C "i^^^ <C |e-B|, it can be shown that for large enough magnetic field and depending 
on the angle 9, a qualitative change occurs in the Coulomb-like potential Vi{R,9); Whereas for = 0, vr 
the potential is repulsive, it exhibits a minimum for angles 9 G]0,7r[ and distances R < 0.005 fm. The 
position of this minimum is proportional to and the depth of the potential at Rmin increases with 

the magnetic field. The exact value of the strong magnetic field will be determined in Sec. V. We interpret 
the appearance of such a minimum as a possibility for bound state formation in D = 4 dimensions. A 
rigorous proof of the bound state formation in the above potentials Vi and V2 is the subject of a separate 
investigation and is beyond the scope of this paper. ^ 

The organization of the paper is as follows. In Sec. II, a brief review of QED in a strong magnetic 
field containing some important results from [4] will be presented. In Sec. Ill, the static potential of 
QED in a strong magnetic field will be calculated perturbatively by determining the vev of the Wilson 
loop of a static fermion-antifermion pair in two regimes of the dynamical mass |q|| <C m'^y^ <^ |ei?| 
and m'^y^ <C |q|| <C \eB\ in the LLL. In Sec. IV a semi-classical Born approximation will be used to 
determine the same particle-antiparticle potentials. In Section V, a qualitative analysis of these potentials 
will be performed and the role played by the angle 9 will be discussed in detail. Sec. VI summarizes our 
results. 



II. QED IN A STRONG MAGNETIC FIELD 



In this section we will briefly review the characteristics of fermions and photons in a constant external 
magnetic field. To this purpose, let us start with the QED Lagrangian density 

C = -\f^uF^^ + ^p-f^' {id^ + eA^) i; - mijiP, (II. 1) 

where the vector field A^ = + where is an Abelian quantum gauge field, F^^, is the cor- 

responding field strength, and A'^^' describes an external electromagnetic field. In this paper we will 
choose the symmetric gauge for A'^^' , i.e., 

^r = f (0'^2,-xi,0). (II.2) 

^ A nonperturbative analysis of the corresponding Schrodinger equation describing the Nambu-Goldstone modes and arising 
from a Bethe-Salpeter equation for bound states shows that at least one bound state can be formed in the attractive 
potential in D = 4 dimensions (for more details, see the explanation in the last paragraph of Sec. II). 
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This leads to a magnetic field in X3 direction. Prom now on, the longitudinal (0, 3) directions will be 
denoted by xy and the transverse directions (1, 2) by x_|_. Using the Schwinger proper time formalism [ol, 
it is possible to derive the fermion and photon propagator in this gauge. As for the fermion propagator, 
it is given by 



SF{x,y) = exj>{'-^{x-yrAf-{x + y)]S{x-y) 



exp 



a,b= 1,2. 



(11.3) 



Here, the first factor containing the external gauge field, A'^ , is the Schwinger line integral [£]. The 
Fourier transform of the translational invariant part S{x — y) reads 

00 



S{k) 







exp I is \_k 



X 



II eBscot{eBs) 

I (^m + 7II • kj|) (1 - tan(eS5)) - 7^ • (l + tan^{eBs)) | , (II.4) 



where ky = (kQjk^) and 7|| = (70,73) and k_L = {ki,k2) and 7_l = (71,72)- After performing the integral 
over s, S{k) can be decomposed as follows 



S{k) =ie"^^(-l) 

n=0 



Dn{eB,k) 



fcy — m? — 2\eB\n^ 

with Dn{eB,k) expressed through the generalized Laguerre polynomials L. 



(11.5) 



Dn{eB,k) = (7II -k,! +m) |2 O 



Ln {2p) - Ln-i {2p) 



+ 47^-A;^Li_i(2p)|. 



(11.6) 



Here, we have introduced p = |^ and 



O = - (1 - i7S'sign(eS)) . 



(11.7) 



Relation pi.SP suggests that in the IR region, with |k|||, |k_|_| <C ^J\eB\, all the higher Landau levels with 
n > 1 decouple and only the lowest Landau level (LLL) with n = is relevant. In strong magnetic field 
limit, the full fermion propagator (jH.SP can therefore be decomposed into two independent transverse 
and longitudinal parts ^ 



SF{x,y) = S'|j(x|j - y||)P(x_L,y±) 



with the longitudinal part 



5'||(x|| -y||) 



„ik||-(x-y)ll 



iO 



7II . k|| — m 



and the transverse part 



\eB\ 



■ exp 



ieB „h n h k-Bl , .9 



a,6 = 1,2. 



(n.8a) 



(n.8b) 



(II.8c) 
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The photon propagator T>^y of QED in an external constant magnetic field in one-loop approximation 
with fermions from lowest Landau level (LLL) is calculated explicitly in 0, 10]. It is given by 

W f,^-S^ + iA + iizZ^Ml!^.,M!L (119) 

where ^ is an arbitrary gauge parameter. Since the LLL fermions couple only to the longitudinal (0, 3) 
components of the photon fields, no polarization effects are present in the transverse (1,2) components 
of T>^y{q). Therefore, the full photon propagator in the LLL approximation is given by the Feynman-like 
covariant propagator Q] 

W^M = ^, (ILIO) 



g2 + q2n ( qi|,qi 



with n(g^,gy) having the form 11] 



^,2 2^ 2a\eB\Nf -^,( \ , 

2 

Here N j is the number of fermion flavors and q = where e is the running coupling. Further, Hi^z) in 
(|ILTT]) is defined by 

g(.).-^=L=ln( ^+^ )-l. (11,12) 

Expanding this expression for ]q|] ^ ''^dyn ^ \^^\ ^'^'^ "^dyn ^ l^yl ^ l^-^l leads to 

„2 



n(qi,qf)^ + ' U for ]q^] « m^^, « ]eS], (ILIS) 

" ^^^dyn. " 

, 9 n, 2a\eB\Nf n , 9, , , , , 

n(ql,qj|)^ L-J-e ^ for m^y„, < ]qj] <C ]eS]. (ILM) 

^ q| 

In [J], it is shown that the kinematic region mostly responsible for generating the fermion mass is with the 
dynamical mass, nidyn., satisfying m^y^ <C ]q|] <C \eB\. Plugging (|II.14p in the full photon propagator 
pi.lOP and assuming that \c[\_\ <^ \eB\, we get 

^-(')^-A' ™* M^ = ^^. (11.15) 

The appearance of a finite photon mass is the result of the dimensional reduction 3 + 1^1 + 1 in 
the presence of a constant magnetic field. This phenomenon can be understood as a reminiscent of the 
Higgs effect in the 1 + 1 dimensional Schwinger model, where the photon acquires also a finite mass. 

As for the dynamically generated fermion mass mdyn., it can be determined by solving the corre- 
sponding Schwinger-Dyson (SD) equation in the rainbow (ladder) approximation, where the effects of 
dynamical fermions are neglected. In this approximation, the gauge invariant dynamical mass is shown 
to have the form 

, / vr / vr \ i/2\ 

rridyn. = CVeBexp 1-- l^—j 1, (11.16) 
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where the constant C is of order one.^ In the improved rainbow approximation, however, the expression 
for nidyn. takes the following form 

,«.,,..=CVliBTFMexp(- ^,^^^';^^^^^ ). (11.17) 

where F{a) ~ (Nja)^/^, Ci ~ 1.82 ± 0.06 and C - 0(1). 

The dynamical mass, rridyn., plays the role of energy in the two-dimensional Schrodinger equation 

(-A + m^j,„. + y(x))*(x) = 0, with A = -^ + ^. (11.18) 



Here, the attractive potential has the form 



y(x) = ^ e.p (S^) Ei (-S^ I . (11.19) 



vr 



and behaves as 

V{.) (j, + In ^) , ^ 0, (11.20) 

where "Je — 0.577 is the Euler constant. To determine the long and short range behavior of the potential, 
the asymptotic behavior of Ei(x) is used [13]. In this context, the problem concerning bound state 
formation is reduced to finding the spectrum of bound states of the Schrodinger equation pi.lSP with 
the attractive potential 01. 19^ . In [4], it is shown that in L> = 4 dimensions at least one bound state can 
be formed. Figure 1 shows the behavior of V'(x) for different B = 10^, 10^, 10^, 10^ and B = 10^^ (from 
left to right) as a function of R.'^ According to this result, whereas for B < 10^ at each given distance 




0.01 0.02 0.03 0.04 0.05 0.06 

R (fm) 



FIG. 1: Potential V{x.) from pi.lOp of a particle-antiparticle pair for different magnetic fields. According to this 
result whereas for B < 10^ at each given distance the absolute value of the potential increases with increasing the 
magnetic field, the shape of the potential does not change for B > 10^. 



^ The problem of gauge invariance of chiral symmetry breaking induced by the magnetic field found in the ladder QED is 
investigated in 12]. 

* Note that discussion here is purely qualitative. It is nevertheless possible to give the strength of the magnetic fields in 
Gauss through the relation 6.8 x lO^^B ~ in GauB [13, H^. Using this relation the potential V{R) is given in GeV=10^ 
eV and the distance R in fm=10~^^ m. 
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the absolute value of the potential increases with increasing the magnetic field, the shape of the potential 
does not change for B > 10^. 

In the next section we will determine the static Coulomb potential between a fermion-antifermion pair 
in two aforementioned regimes in the LLL approximation by calculating the vev of the Wilson loop 
perturbatively. Note that although our method is a perturbative one, our results depend explicitly on 
'n^dyn. which is determined nonperturbatively in the literature 0|. 



III. THE WILSON LOOP AND THE MODIFIED COULOMB POTENTIAL IN A STRONG 

MAGNETIC FIELD 



In ordinary Quantum Field Theory with no background magnetic field, the Wilson loop appears as one 
of the most efficient tools for probing the large distance properties of strong coupling QCD. It provides 
a natural criterion for confinement through the area law. In contrast to QED, where the field lines 
connecting a pair of opposite charges are allowed to spread, one expects that in QCD the quarks within 
a hadron are the sources of a chromoelectric flux which is concentrated within narrow tubes connecting 
the constituents. Since the energy is not allowed to spread, the potential of a quark-antiquark pair (qq) 
will increase with their separation, as long as vacuum polarization effect do not screen their color charge. 
This picture of confinement can be checked by computing the nonperturbative potential between a static 
quark-antiquark pair in the path integral formalism. 

It is the main purpose of this section to investigate the properties of the Wilson loop concerning the 
bound state problem of QED in the presence of a strong magnetic field in the LLL approximation. Here, 
the magnetic field plays the role of a strong catalyst and even the weakest attractive potential between 
fermions is enough for dynamical mass generation and bound state formation. Before calculating the 
modified Coulomb potential by determining the vev of the Wilson line of a static fermion-antifermion 
pair in two regimes of dynamical mass, |qy| <C m'^y^ <^ \eB\ and m'^y^ <C |qy| <C \eB\, we briefly 
review the path integral formulation leading to the Coulomb potential between slowly moving particles 



in ordinary QED. Here, we keep closely on the notations of 16(] and the classical review article [a]. 

Imagine creating a fermion-antifermion pair, at space time point = and adiabatically separate 
them to a relative distance R. This configuration will be hold for a time T — > oo. Finally we bring t/jip pair 
back together and let them annihilate. The Euclidean amplitude for this process is the matrix element 
of the Hamiltonian evolution operator e~^'^ between the initial and final states, \i) and |/), respectively 

(z|e-^^|/). (IILl) 

Here, \i) and |/) represents a ipip pair a distance R apart and H is the Hamiltonian (for more details 
concerning the exact mathematical structure of the initial and final states, \i) and |/), see [igI). In the 
path integral representation the matrix element pil.ip can be expressed by 

where we have already skipped the subscript E for the Euclidean action S. Here, the current should 
describe the closed worldline C of the creation and annihilation of ^ipip pair. For a closed worldline of the 
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heavy fermions ()III.2p simplifies to [l6| 



— S+ie ^ Afj,{x)dxn 

Since in a closed path \i) and |/) are identical and since the process is static, the Hamiltonian H is purely 
potential and the left hand side (l.h.s) of pil.sp is 

(i|e-^^|/) = e-^(^)^. (IIL4) 

Here, V{R) is the fermion-antifermion potential. Taking the logarithm of (|III.3P by plugging pil.4p on 
its l.h.s, it is given by 

V{R) = -l\m lln(T^c[^]), (111.5) 



where we have introduced the "loop-correlation" function - the Wilson loop 



Wc[A]=e c , (III.6) 



and used the definition 



2 

For S being the QED action, the Coulomb potential V{R) = —-^^^ can be analytically calculated in 
the quenched approximation, i.e. when the vacuum polarization effects arising from the presence of 
dynamical fermions are neglected. To calculate this Coulomb potential, one expands (VFci^]) in powers 
of the background field , to get 

{Wc[A]) = {1 + ie £dx^A^{x) - ^ £ £dx^dyuAf,{x)A^{y) + ■ ■ ■) . (III.8) 

Here, the second term including only one gauge field does not contribute. As for the other terms, the 
terms with an odd number of external photon lines does not contribute to the above expansion. This 
is because of the Furry's theorem, that holds in ordinary QED in contrary to QED in the presence of 
external magnetic field ^. Plugging pH.Sp in pn.5p . the Coulomb potential is given by 

V{R) = lun^ In ^1 - y jf jf dx^dy^ D^^{x, y) + ■ ■ 

= hm ^ / / dx^dy^ D^^{x,y) + 0{e^), {111.9) 



where we have expanded the ln(l — • • • ) for weak coupling constant e. The integrand is the photon 
propagator in the coordinate space 

D^,{x, y) ^ {A^{x)A,{y)) = (HI.IO) 



^ The phenomenon of photon sphtting 7 ^ 27 in the presence of external magnetic field [2Q| is the best counterexample for 
the validity of Furry's theorem for QED in external magnetic field. 
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To compute the double integral (|III.9p . one chooses a rectangular contour. Because, in the Euclidean 
space, the integrand pil.lOp is proportional to 6^,^, the double integral receives a contribution when x 
and y are located on the segments of integration contour which are parallel to each other. The double 
integral therefore reduces to 



V{R) = lim — / dTi / dT2 DooiR,Ti-T2 



o2 rT [-T 



JO 



lim — / dTi dT2 ^. (III. 11 



47r2[i?2 + (Ti-T2)^ 



Here R= |x — y|. Performing now the double integration, one arrives first at [. 

and then after taking the limit T — > oo, at the Coulomb potential of a fermion-antifermion pair 

V{R) = (III.13) 

After this brief excursion, let us now turn back to our original problem of computing the modified 
Coulomb potential of QED in the presence of external strong magnetic field. In the regime of lowest 
Landau level (LLL), the full photon propagator D^i, in the momentum space is given by (jll.lOp . where 
n(q|, q^) is defined in pi. lip . It is indeed a very difficult task to calculate the potential of static fermion- 
antifermion pair for full photon propagator pi.iop . We will therefore consider only two different regimes 
of dynamical mass |qy| <C <C |ei?| and m'^yn |q|| |ei?|, with n(qy,q^) given in pi.isp and 

pi.l4p . respectively, and will calculate the modified Coulomb potential V{x) in these two regimes. 



A. Modified Coulomb potential in |q|| ^ nT-dyn ^ l^^l regime 
Substituting (flLTSl) with Nj = lm (ULToI) the full photon propagator in this regimes reads 

To calculate the Coulomb potential by pil.op . we need the LLL photon propagator in the coordinate 
space, i.e. the Fourier transformed of piI.14p 

~ II f d'^q e*^^' 

'D.Ax) = -^9l y (2^)4 2 r^T^' ^"^'^^^ 
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After a lengthy but straightforward calculation of the above integral over q (see Appendix A for more 
details), we arrive at the LLL photon propagator in the coordinate space 



6, 



7 4a2 + 7ii^ sin^ 6 802^?^ sin^ 9 



l + - ^-TTT^ + 



ai 2Pai ' 2p^al 



+ ^ 1 ^r. + 



2/3ai 8/32af 
' ' ' 2K sm & -\ H ——7 it sm + ' 



(III.16) 



\eB\(3aj \ (3ai \ 2 j J p^aj \ 16 1 J S-fp^af 

Here, 



are constant c-numbers, and Oj = 9,T) , i = 1,2 are defined by 

ai{R,9,T) = + R^f{a,9), with /^(a, 0) = 1 + sin^ 0, 

a2{R,9,T) = /?7 (r^ + i?2 cos^ 0) . (III.18) 

In all the above expressions R = |x| is the distance between the static fermion and antifermion pair, 
is the angle between the particle-antiparticle axis and the direction of the magnetic field B, i.e. the X3 
direction, and T = ixq is the Euclidean time. Note that here, to determine the LLL photon propagator, 
we have to use of the approximation <C \eB\, which is valid in the regime of strong magnetic field. 

To determine the static potential between the Tpip pair, we have to calculate the double integral [see 
pil.lip leading to the expression of the ordinary Coulomb potential] 

V{R,9)= lim r dTi r dT2Voo{R,9,Ti-T2), (11119) 

which can be easily simplified to 

V{R,9) = -2e^ I dT Voo{R,9,T). (IIL20) 
Jo 

Using now the definition of T>fj_u{R, 9, T) in piI.16p . and performing the integration over T, the modified 
Coulomb potential in the regime |qy| <C 'm'^yn can be given by 

m,«) = -|(^.(a.<')-^^^ + ^^^), (111,21) 
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with 



1 

T 

1 Zcos^e ?>sm^e 15 sin^ 61 005^6' 



, , 9 / Ssin^ei 35sin^6l 

^3(a, 0) = +— ^ 1 - — ^ + 



16/5 V 4/3/2 128/32 /4 
Ibcos^e 7 o„ „ . 9 63sin2(9/^ 3sin2(9\ 693 cos^ 6* sin"^ 6* 



8/7 V" 4/3/2^"''""" " ' ' 8/32/4 y^--™ — iQ J 2 56/32 /6 

3/5 ■ 35sm20 / sm2e\ 315cos20sm4| 

(111.22) 



In Sec. V we will present a qualitative analysis of the above potential emphasizing the role played by the 
angle in a possible bound state formation. 



B. Modified Coulomb potential in ^ |q|| ^ \eB\ regime 



To compute the inter-particle potential in this regime, we substitute (IILT4|1 in (flLTOl) . For Nj = 1, the 
full photon propagator in the coordinate space reads 



(2^)%2_«exp 



'2IZBI 



(III.23) 



The integration over q can be performed using the approximation <C |ei?| which is valid in the regime 
of LLL dominance. After a straightforward calculation (see Appendix B for more details), the propagator 
is given by 



where 



C 



-'fJ.U 



c 



47r2 (1 - f ) ^/WTlWe)' 



2a\eB\ 



vr 



and g'^{e) = cos^e + 



sm 



(111.24) 



(111.25) 



The fermion-antifermion potential in this regime is then calculated using the expression ()III.19P or equiv- 
alently (IIII.20p and reads 



To evaluate the integral we make use of 

7-2^+1 



dT 



(III.26) 



dT- 



Choosing t = T , z = Rg{0), u = 1 and /i = — ^, we arrive at 

v{R,e) = - 



^^;^^K._,_i(Cz), C > 1, Re(^) > -1. (III.27) 



a 



{l-^)g{e)R 



(III.28) 



11 



^ and r(^) = -y/vr. Apart from its 6 dependence our result is 



Here we have used Ki ix) = 

2 

comparable with p.ip from [7J. 

The potential ()III.28P is indeed comparable with the ordinary attractive Yukawa potential 



ct 

Yukawa (-^) = ^6 



-mR 



(III.29) 



where a 



and the effective photon mass m — > M^ff {9) = C,g{0). This result is also in agreement 



with the general result about the massive photons in a strong magnetic fields in the LLL approximation. 
As we have seen in Sec. II, in the regime rn^y^^ <C |q^| ^ the 3 + 1 dimensional QED in the LLL 



approximation is reduced to a 1 + 1 dimensional Schwinger model where the photon acquires a finite 



mass = y/ ^"!*^^!^/ ^ j\r^ the number of fiavors [see (jll.lSP for more details]. For Nj = 1, we have 
My = Thus the effective photon mass in this regime is given by M<,ff, = M^g{6) and depends explicitly 
on the angle between the particle axis and the direction of the external magnetic field. Figure 2 shows 
the ^-dependence of g{6) which is understood to be the ratio For all values oiO £ [0, vr] the effective 

mass Mcfl and the aforementioned photon mass are almost equal. Here a is fixed to be a = 1/137. A 



qualitative analysis of the potential ()III.28p emphasizing the role played by the angle 9 will be presented 
in Sec. V. 



4eff/MY 




0.5 1 1.5 



2.5 3 



FIG. 2: The ratio ^ - 



g{6) for 6 G [0,7r] in the regime m'^ ^ \qf, \ ^ \eB\ of LLL dominance. 



In Sec. IV we will use the Born approximation to determine the static Coulomb potential in two regimes 
of dynamical mass in the LLL approximation. Eventually we will study the behavior of the potentials 
pil.2ip and pil.26p as a function of the angle 9 for two limits of large and small distances. 



IV. BORN APPROXIMATION AND THE STATIC COULOMB POTENTIAL IN A STRONG 

MAGNETIC FIELD 

In the nonrelativistic quantum mechanics the relation between the scattering amplitude A4 and the 
potential is given by the Born approximation 

{p'\iM\p) = -iV{ci){27T)6{Ep, - Ep), (IV.l) 

where p (Ep) and p' [Epi) are the momenta (energy) of the incoming and outgoing particles, respectively, 
and q = p' — P- For ordinary QED with no background magnetic field, for instance, the amplitude of a 
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particle-antiparticle scattering is given by 2l| 



iM ~ - , (IV.2) 

|p-p'r 



Comparing with (|IV.ip . the attractive (classical) Coulomb potential V{q) is thus given by 

^(q) = with |q| = |p - p'|. (IV.3) 

After a Fourier transformation into the coordinate space, the same potential reads 

where R = |x|. Furthermore, to include the quantum correction into the result, the modified Coulomb 
potential can be calculated from 

'-W--/|?-" %Mi-n(,^)) - 

where n(q) in the ordinary QED is defined by the vacuum polarization tensor 

n^.(g) = ('zV'^-'zV)n(g2), (iv.6) 

and is given by 

— / dx x{l-x)log(— ^ 

Choosing qo = and plugging this relation into pV.5p . after some straightforward calculation [21,], one 
arrives at the so called Uehling potential 

-2mR 

r 



n(9') = / dxxil-x) log — ^ . (IV.7) 



In this section, we will use the above method to determine the potential of QED in a strong magnetic 
field in two regimes of dynamical mass |q|| ^ "i^^^ <C \eB\ and ^ |q|| <C |ei?| separately. We 

will show that our results coincide with pil.2ip - pil.22p in the regime |qy| <^ ^dyn ^ \^^\ ^^'^ with 
pil.26p in the regime <C |q|| <^ \eB\. 

A. Modified Coulomb potential in |q|| ^ ?t7,^^„ <C \eB\ regime 

We will start using the relation pV.5P where n(q) in this regime is given by pi.isp . The modified 
Coulomb potential in this regime is therefore given by 

q2 + 2^ g| e 2ks| 



where we set qo = 0, and the factor 7 = "j^^^ is defined in piI.lTp . The 3-dimensional integral over 

dyn. 

q can be performed by making use of the same methods as was shown in Appendix A. The main steps 
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of the calculations are as follows. First using the relation ()A.2p and (jA.Sp . we write the integrand in the 
form 



y(x) 



ds 



dq±q±dq3dip' i(q^Jicos e+q^Rsmecos{<fi-ip')) 
(2^)3 



exp -s + 93 + 



(IV.IO) 



where 6 is the angle between the particle-antiparticle axis and the external magnetic field. The integration 
over If' an can be performed using (fAlHI) and (lAl2]) . respectively. We arrive therefore at 



ds\ / — 



vr f dq±q± 



2 R'' cos^ e 



Jo {q^RsmO) e-'^^^e'-^ 



(IV.ll) 



where we have introduced 



7e5 „ 1 
-e 



(IV.12) 



to simplify our notations. At this stage, to perform the integration over and eventually over s, one 
can use the approximation |q^| <C |e-B|, which is valid in the LLL approximation. The factors ^ and 
are therefore given by 



4^2 8|e5|^ 
37^ 



qi 



7 



1 / 2 / 7 2 



ql 



where ^ = (2 + 7|ei3|). Using these expressions, the potential is given by 



y(x) 



X U + 



ds\ 



Tqi 

4^ 



2iT f dq±q± 



J (27r)2 
^1 



g-'^iie Jo {q±R sine) 



si 



+ ql 



72 f3 36i^ 62 



16^2 ^2 2s2^2 y I6^|e5 



7 



(IV. 13) 



(IV. 14) 



where 61 = R^cos^O. The result is similar to ()A.16p and indeed we have used the same approximation 
which was done in (lAlSj) and led to (|AT6)) . Using now the integrals li, i = 1,2,3 from ()A.2ip - (|A.23p . 
pV.14p can be written as 



F(x) 



ds I2tt _i 
e 



Tqi 

4^ 



(27r)2 V si 
h 



1 



si 



7' f3 36i ^ 6? 



16^2 ^2 si 2s2^2 y iQ^\eB 



7 



Here, the integration over s can be performed using (|A.25p . After replacing 

2i?2 cos2 e 



R^ sin^ e + 



APR^f{9), 



i 



(IV.15) 



(IV. 16) 
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where /3 = (2^) ^ and f'^{a,6) are defined in piI.lTp and (jlll.isp . we arrive at the same potential 
pil.2ip - pil.22p . The potential has therefore the general form 



ViR, 9) = -- [A,{a, 0) - + ^4 J ' 

with A, i = 1, 2, 3 from pil.22p . 



B. Modified Coulomb potential in rri^^y^ ^ |q|| ^ |ei3| regime 

Starting from PV.SP and plugging n(g) corresponding to the relevant regime rn^y^ <^ |q|| <^ \eB\ from 
pi.l4p and choosing qq = 0, we arrive first at 

— ^. (IV.IT) 



where C = y — ^f— ^ is already defined in pil.25p . To perform the three dimensional integral over q, 
we will use the same methods described in Appendix B, to evaluate the full photon propagator in the 
coordinate space. First using Schwinger's parametrization technique ()B.2p . the potential reads 

F(x) = -e'Jdsj '^'^^l^^^^f'^' e-^(''^^«-^-^(^-^')+'?3«™^^) exp (^-s + ^ e'^j j , (IV.18) 

where we have used (jA.Sp to bring q • x in a useful form. Now integrating over ip' and using (jA.lOp 
and (|A.12p . we arrive at 

oo 

F(x) = -e2 Jds^J^ e-^Jo(gxi?sin^)e-<%-(i"fK, (IV.19) 





|2 



where we have used the IR approximation |qp ^ |ei?| relevant in the LLL regime 



Next, using (jB.lOp the integration over q± can be performed. The potential is therefore given by 



oo 



y(x) = -j^ — — / — e ~ , with g^{9) = cos^ 6 + 







The ^-dependent function g(6) is defined already in pil.25p . Finally defining a new variable s' = C^s, 
and using (1B.13P the potential F(x) in the riT^y,^ <^ |q|| ^ \eB\ in the LLL dominance can be calculated 
and reads 

V(R 0) = - e-^3im 

which is the same potential as pil.28p which was found using the Wilson-loop technique. 
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V. A QUALITATIVE ANALYSIS OF MODIFIED COULOMB POTENTIALS IN THE LLL 



In this section we will consider the potential ()III.22|1 and pil.28p 



ViiR,i 



a 
R 



Al[a,B) h 



i?2 



V2{R,9) 



a 



(1 - ^)9{0)R 



in the first |qy | ^ ''^dyn ^ k-^l) second m'^y^ <^ jq^ | <^ \eB\ LLL regimes, respectively. We will 

eventually compare these potentials with the modified Coulomb potential p.ip . 

As for Vi{R,9), Figure 3 shows this potential for different choices of the magnetic field B = 
10^ 10^ 10^, 10^ (from right to left) and different 9 = 0, vr/S, 2tt/3, it. 



Potential Vi (R,^?) in qn << md << I eB | Regime for different 6 and B 
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FIG. 3: Potential Vi{R,d) for different B and 6. No qualitative changes occurs by varying the angle 9. 



According to this result, for i? — > the potential falls more rapidly to —oo the smaller the magnetic field 
is. Although the scales in which the potential Vi is plotted for 9 = 0,it are different from the scales in 
which it is plotted for 9 = 7r/3,27r/3, but their shapes are almost the same, i.e. no qualitative changes 
occurs by varying the angle 9. This situation changes by neglecting the coefficient 7^.43 comparing 
with Ai and 7^2 in yi{R,9). Figure 4 shows the behavior of the coefficients ^1, 7^2 and j'^As of the 
potential Vi as functions of the angle 9 for different mag netic fields B = 10^10^10■^ and B = . 
Choosing a = 1/137, it turns out that the coefficients ^1,7^2 and 7^.43 are positive V0 € [0, vr] and for 
any choice of constant magnetic field B. However, as it is shown in Figure 3, .A3 decreases rapidly with 
increasing magnetic field. Whereas for B = 10^ two coefficients ^.1,7^.2 and 'y'^A3 are comparable, for 
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The Coefficients A1/10-^ yAj/IO"' and /Aa/IO"^ for B=10^ 
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The Coefficients Ai/10"', yAs/IO"'" and /Aa/IO"^^ for B=10' 
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The Coefficients A1/10-^ rA2/10"'^ and /Aj/IO"'^ for B=10^ 
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The Coefficients Ai/10-'',yA2/10"''' and /Ag/IO"^" for B=10^ 
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FIG. 4: Coefficients Ai, 7^2 and j'^As for different magnetic fields. As it turns out, ^3 decreases rapidly with 
increasing magnetic field. 



B = 10^, j^As is 10^ times smaller than ^1, and for B = 10^ this difference is ~ 13 order of magnitude. 
We conclude therefore that for strong magnetic field B > 10^, the coefficient 7^.43 in Vi{R, 9) is negligible 
comparing with Ai and 7^2- Thus for B > 10^ the potential Vi{R,e) can be replaced by 

F3(iJ,#)=-^(.4,(«.«)-^«), (V.l) 

which has its minimum at 



Rmin{B,9) — 




(V.2) 



Figure 5 shows V3(i?, ^) for different choices of the magnetic field B = 10^, 10^, 10^ and 10^ (from right 
to left) and different = 0, 7r/3, 7r/2, 27r/3, tt. Whereas for = 0, tt the potential is repulsive, it exhibits a 
minimum for angles 9 G]0,7r[ and distances R < 0.005 fm. The depth of the potential at Rmin increases 
with the magnetic field. We interpret this effect as a possibility for boTind state formation. The answer 
to the question concerning the existence and the number of bound states in the potential V3 is beyond 
the scope of this paper. 
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FIG. 5: Potential V3{R, 6) for different B and 6. Bound states can be formed for 9 e]0, 7r[ and for strong magnetic 
fields B > 10^ in the regime R < 0.005 fm. The depth of the potential at Rmin increases with the magnetic field. 



In Figure 6 the behavior of Rmin from (|V.2p for different 6 is studied. As it turns out, for different 6, the 
position of the minimum of the potential is proportional to \j\fB. 




lo'sBslo' 



FIG. 6: Behavior of the minimum of the potential V3 as a function of 10^ < B < 10^^ for different < 9 < tt/2. 
For 7r/2 < 9 < t:, Rmin shows a symmetry in changing 9^9 + tt/2 (see the 3-dimensional figure on the r.h.s.). 
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Potential V2(R,fl) in m/« « | eB | Regime for different 8 




FIG. 7: Potential V2{R,0) for different B and 6. No qualitative changes occurs by varying the angle 9. This 
potential is comparable with the potential (jl.ip from 7] in the scaling regime. 



Let us now consider the potential V2{R,6) in the second regime m^y^ <C |qy| <C \eB\ of LLL. Figure 
7 shows its behavior for different magnetic field B and angle 9. Again no qualitative changes occurs by 
varying the angle 6. As it is pointed out in the introduction, this Yukawa-like potential is comparable 
with the potential (jl.ip from [7|] in the scaling regime. 



VI. SUMMARY 



In this paper the static potential of QED is calculated in the presence of a strong but constant magnetic 
field using two different methods. First a perturbative Wilson loop calculation is performed for two 
different regimes of dynamical mass |q|| <C rn'^y^ <C \eB\, and m^^„ <C |q|| ^ \eB\ in the lowest 
Landau level (LLL). The resulting potential is then compared with the potential arising from a modified 
Born approximation. The results coincide. Comparing with the recently calculated potential of point like 
charges in 0], our potential shows a novel dependence on the angle 6 between the particle-antiparticle axis 
and the direction of the magnetic field. A qualitative analysis of these modified potentials is performed in 
the previous Sec. V. As for the potential pil.22p from the first regime |q|| <C rri^y^ <C \eB\, it turns out 
that for strong magnetic field B > 10^, the coefficient of the term, i.e. j'^As, is negligible comparing 
to the both other coefficients .Ai and 7.A2 from the expected Coulombian and the additional 
terms. Keeping only these two coefficients, a qualitative change occurs in the Coulomb-like potential 
which depends on the angle 9; Whereas for 9 = 0,11 the potential is repulsive, it exhibits a minimum 
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for angles 6 G]0,7r[ and distances R < 0.005 fm. The depth of the potential at Rmin increases with the 
magnetic field and the position of the minimum is proportional to We interpret this effect as a 

possibility for bound state formation. A rigorous proof of the existence and the number of bound states 
in the above potential is the subject of future investigations. 

As for the potential pil.28p of the second regime m^^^ <^ |q|| <C \eB\, it is comparable with a screened 
Yukawa potential. It is well-known that in this regime the photon acquires a finite mass proportional 
to y^|e-B|. According to our result, the photon receives also a modified effective mass depending on the 
angle 6 between the particle axis and the direction of the magnetic field. But, this dependence is indeed 
negligible for fixed a = 1/137. This result is in good agreement with the result recently found in [3]. 

It would be interesting to study the renormalization group (RG) improvement of these potentials, as 
they both depend on the coupling constant e, which is taken to be a bare parameter in our calculations. 
Using the RG improved potential it is possible, according to [8], to determine the Callan-Symanzik 
/3-function of QED in the presence of strong magnetic field in the LLL approximation. 
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APPENDIX A: THE LLL PHOTON PROPAGATOR IN |q|| < mly^ < \eB\ REGIME 

In this section we will perform the integration over q in (|III.15jl to deter mine the full LLL photon 
propagator in m^y^ <^ |qy| <C \eB\ regime (IIII.16h . To start, let us give the integral pil.lSh in the 
Euclidean space 

where 7(a) = ^" — . Here, the Euclidean coordinates xq = —ix4 as well as qo = —iq^. We have 
further used the notations q = ((71,^2,93) and x = (xi, 2:2, 2:3). The scalar product is defined therefore 
by q • X = Yl^i=i Qi^i- Using the Schwinger parametrization technique 



00 







ds e-"' = -, (A.2) 
a 



the above expression (jA.ip can be given by 

foo r / / ^Y^Ai^Di / q2^ 

'2kB| 



^..(-) = Si [ ds J ^e-(--+-) exp (qj + q V (,i + qj) exp (-^ 



(A.3) 



Let us first evaluate the integral over ^4 



+00 



^^46"^"*^* exp ( -sql { 1 + ] ] ^ ^^ ^^ 
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For the new variable 

s' = s 1 + ^ , (A.5) 

we get 

^^^g-i54X4-s'5| ^ J-e-i7. (A.6) 

'OO V 'S 

To integrate the variable ^3 component, let us go into the polar coordinate system. Here, the product 
q • X is given by 

q • X = qR sin 9 sin 6' cos ((/? — (/?') + qsR cos 0' cos 9, (A. 7) 

with g = |q| and R= |x|, (9') the angle between x (q) and the external magnetic field, which is assumed 
to be in the xs-direction. Now using ^3 = qcos9' and q_L = qsm9', ()A.6P is given by 

q • X = q±Rsm 9 cos (99 — ip') + q^R cos 9. (A. 8) 

In these cylindric coordinates the photon propagator (jA.Sh is given by 

2 \ \ \ .2 



2,2, 7(a)|e5| 2 \\\ 



xexpi^-si^qi + cii + ^ 1^3^ exp (^-^ji^J J J ^ -e 4.', (A.9) 

where we have inserted the expression from ()A.6j) . Now using the integral representation of the Bessel 
function Jq 

f-2TT 

27r 

and performing the integration over if we get 



Jo(z) = ^ I dip'e-''^°<^'-^\ (A.IO) 



^..(x) = Si ds I !^^gi^e--Hcos.-.',3Vo iq^Rsin9) e'^^l^e-i. (A.ll) 
Performing now the integration over 53 using 

" d<^3e-^««-^^-Ai = (A.12) 



The photon propagator (jA.lip therefore reads 



~ II ; r/s fr/flifli 9 (izi|+fl2 cos2 s) 

= 5l 4 / ^^Jo(g±i?sm0)e-^'iie-' (A.13) 



To perform the integration over q± we proceed as follows. Using first the IR approximation qj^ <^ |ei?| 

q2 , „4 

in the regime of LLL dominance, and expanding the factor e ^ksi ~ 1 — -|- the relation 

s'/s from (jA.Sp can be written as 

^ = ^fl + ^+fl4-^^lqiV ^ith e^(2 + 7|ei3|). (A.14) 
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Plugging this result in (|A.13P the photon propagator is given by 







J 47r2 "^^^ > \ 2i ^i\eB\ 



Again using the IR approximation and expanding the exponent we get 
= / / dg^gx Jo (g±i2sin0) e"^'!^ 

471" Jo ■S JO 



( ^ ^1 \ , „4 



where 

„2 , d2 „„„2 



61 = X4 + i? cos 
To perform the integration over we use [3] 

Jo pk'^t{u + i) 2 '2 V^^y 

where Mm (y) is the Whittaker function defined by 



ME^^{y)=y 2 e 2$( ,u + l;y] 



where 



'^^ ~ ^ rl! ^ r(T + l) 2! ^ . 
For the first term on the second line of ()A.16P we get therefore 

h{s)= / dq±q_L Jo {q±R sin 6) e '"^^ = —e 4s , 
Jo 2s 

and for the second term we get 

i?2 sin^ e 



4s 



'igxgi Jo ((Z±i?sin0)e-^'ii = — e'^H^^ h 

whereas the third term yields 

/3(s) = / Jo iq±Rsm9) e = e 4. 1 + 

JO V 

Inserting ()A.2ip . (|A.22p . and ()A.23P in ()A.16P we arrive first at 



32^2 



4^2 \^ si 8s2^2 y 8^|g^| \^ 2si 



his) 



22 



Finally, the integration over s can be evaluated using [l3| 



ds 



r(n- 1) 



Choosing the notation 5o = sin^ 9 + , we get for the first term in (|A.24p 



(A.25) 



for the second term 



ds 







and the third term 



e 4s 



ds _^ 



2^e&o' 



^Ivl ( 46i _^ 126ii?2sin2 6i 1 2i?2 sin^ 6* 



JjJLU 



ds 



4vr2 7o s 



4^2 



h\ 



si 852^2 7 8e|eS 



2si 



his) 



6i?2 sin2 6i?4 sin^ 6 
1 ^ + 



bo 



/)2 



7r2e464 

44^/7 
"7r2|eS|^253 



8 /^„2 • 2^ ftiA 20 /iJ^sin^' 
bo \ hi \ 2 



+ 



1 



^(^i?2gin2 0+|]+|fi?Sin4 + 



861 i?2 sin2 6l 



Cb'o 
60biR^ sin^ 



(A.26) 



(A.27) 



(A.28) 



Further, to find the LLL photon propagator in the regime |q2| <^ ^dyn ^ \^^\ from OILlGh . we will 
first replace 

a2{R,9,X4) 



bo 4/3oi(i?,6',X4) 



bi 



Hi 



and 



(A.29) 



where /3 and a,,? = 1,2 are defined in pil.isp . Adding then the results from ()A.26p . ()A.27p and ()A.28P 
together we arrive at the propagator piI.16j) . 



APPENDIX B: THE LLL PHOTON PROPAGATOR IN m^^^ < |q|| < |eB| REGIME 

In this section we will perform the integration over q in pil.23p to determine the full LLL photon 
propagator in m^^^ <C |qy| <C |ei?| regime (jIIL24p . To start, let us give the integral piL23p in the 
Euclidean space 

where we have introduced the Euclidean coordinates xo = — ix4 as well as qo = —iqA- Here, q = (gi, (?2, 93) 
and X = (xi, X2, X3) and q • x = X]i=i Qi^i- Using the Schwinger parametrization 

ds e-"' = -, (B.2) 
a 



23 



the above expression ()B.ip can be given by 



V,.{x) =6ll ds I e-(^-'^+'^-) exp ( [ qj + + ^^e"^ ] ] . (B.3) 



Performing the integration over (74 



+°° 2 ■ /tt -4 



00 



we arrive at 

P,„W=*;„y^"d.yf.-i/^0e--<.«exp('-Jq^ + ?;^.-^*.)V (B,5) 

To integrate the q component, we follow the same steps as in the previous section [see (jA.7p and (|A.8p ] 
leading from ()A.3p to ()A.9p . In the cylindric coordinates the photon propagator (|B.5p is given by 



sin 9 cos{ip~ ip' )+q'iR cos I 



2 , 2a|e5| 



2 



xexp^-s^q^H ^ -e j j . (B.6) 

Here if' G [0, 27r] and q± = \(l±\ & [0, 00). Now using the integral representation of the Bessel function Jq 
from (jA.lOp and performing the integration over 99 we get 

xexp f-s (qi+(?i + ^^^e~^ I I , (B.7) 



TT 



where we have written q^ = qj^ + (/I- Performing now the integration over in the same way as the 
integration over q^ [see ()B.4p ] we get first 



and then 



with = H£!i£M Xo perform the integration over q_L, we use the approximation, q^ <^ which is 

valid in the regime of LLL dominance. After expanding the exponent 



the integration over q_|_ can be written as 



^-'''l^dq, q, Joiq,Rsine)e~<^~^hl = ^T^I^ exp ("if^) • (B-IO) 
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To evaluate the q± integration, we have used [l^ 

J dz z''^^Jy{l3z)e-^' = (2^)1^+1 ^"^' ^^^^^ > > 0' (^-^1) 

by choosing 7 = s(l — ^),/5 = i?sin0, and = 0. Plugging this result in (IB.9|) we arrive at 



167r2(l-a)i„ »2 

where 5(0) is defined in (|III.25p . Defining a new variable s' = C^s, the s'-integration can now be performed 
using 

„ ;St'-'-- = (!) 2A'.(.)^ (B.13, 
We arrive finally at the full LLL photon propagator in the regime fn^y^ <C |q^| <C \eB\ is then given by 



4-' c 



For the notation X4 ^ T we arrive therefore at our results from ()III.24p . 
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